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Abstract 

This paper describes a novel scheme for integrating a coherent THz source in 
the baseline of the LCLS facility. Any method relying on the spent electron beam 
downstream of the baseline undulator should provide a way of transporting the 
radiation up to the experimental floor. Here we propose to use the dump area access 
maze. In this way the THz output must propagate with limited size at least for one 
hundred meters in a maze, following many turns, to reach the near experimental 
hall. The use of a standard, discrete, open beam-waveguide formed by periodic 
reflectors, that is a mirror guide, would lead to unacceptable size of the system. 
To avoid these problems, in this paper we propose an alternative approach based 
on periodically spaced metallic screens with holes. This quasi-optical transmission 
line is referred to as an iris line. We present complete calculations for the iris line 
using both analytical and numerical methods, which we find in good agreement. 
We present a design of a THz edge radiation source based on the use of an iris line. 
The proposed setup takes almost no cost nor time to be implemented at the LCLS 
baseline, and can be used at other facilities as well. The edge radiation source is 
limited in maximally achievable field strength at the sample. An extension based 
on the use of an undulator in the presence of the iris line, which is feasible at the 
LCLS energies, is proposed as a possible upgrade of the baseline THz source. 



1 Introduction 

The accelerator complex at the LCLS produces ultra-short electron bunches 
approaching sub-hundred fs duration. It is natural to take advantage of 
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these ultra-short bunches in order to provide coherent THz radiation [1]. 
In fact, intense, coherent THz radiation pulses can be produced from the 
sub-hundred fs electron bunches at wavelength longer than, or comparable 
with the bunch length, leading to radiated energy levels proportional to 
the square of the electron number, in contrast to the incoherent case when 
energy in the radiation pulse scales linearly with the number of electrons 
involved in the process. The result is an enhancement in radiation intensity 
of up to 9 - 10 orders of magnitude. 

The exploitation of such kind of coherent THz source as a part of the LCLS-II 
baseline user facility has been proposed in the LCLS-II Conceptual design 
report [1]. THz radiation pulses can be generated by the spent electron beam 
downstream of the X-ray undulator. In this way, intrinsic synchronization 
with the X-ray pulses can be achieved. A first, natural application of this 
kind of the photon beams is for pump-probe experiments. Through the 
combination of THz pump and X-ray probe, LCLS-II would offer unique 
opportunities for studies of ultrafast surface chemistry and catalysis [1]. 
Also, the LCLS team started an R&cD project on THz radiation from the 
spent electron beam downstream of the undulator [2]. The THz is generated 
by inserting a thin Be foil into the electron beam [3]. In this paper we 
describe a novel scheme for integrating such kind of source in the baseline 
of the LCLS facility. 

The transport of the THz radiation from the LCLS beam dump area to 
the near experimental hall constitutes a challenge. A major constraint for 
the LCLS baseline case is constituted by the upstream shielding wall and 
by the Far End Enclosure (FEE) downbeam shielding wall. It follows that 
transport of the THz beam can only be achieved using the dump area access 
maze, and relying on a limited size [4, 5]. However, in this way, the THz 
output must propagate at least 100 meters to reach the near experimental 
hall. Since THz beams are prone to significant diffraction, a suitable beam 
transport system must be provided to guide the beam along large distances 
maintaining it, at the same time, within a reasonable size. Moreover, the 
THz beamline should be designed to obtain a large transmission efficiency 
for radiation over a wide wavelength range. Usually, the focusing of the 
THz beam can only be achieved with reflective optics because lenses made 
from any material would reflect and absorb all radiation. In order to cope 
with the unacceptable size increase, in this paper we propose an alternative 
solution to a mirror guide, based on the use of periodically spaced metallic 
screens with holes. This quasi-optical transmission line is referred to as an 
iris line, or an iris beam waveguide. The eigenmodes of the iris line have 
been calculated numerically for the first time by Fox and Li [6] and later 
obtained analytically by Vainstein [7, 8]. When the Fresnel number of the iris 
line is large, eigenmodes are characterized by rather small diffraction losses. 
For instance, at a wavelength A ~ 0.1 mm, for an iris radius a ~ 5 cm, and 
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for a the distance between the iris b ~ 30 cm, the Fresnel number is given by 
a 2 /(Ab) ~ 10 2 , and diffraction losses of the principal eigenmode are found to 
be about 10% in 100 meters. Also, the iris line is rather stable with respect 
to screens misalignments. When the screens are adjusted in transverse and 
longitudinal directions with an accuracy better than 1 mm, misalignment 
do not result in extra diffraction losses within the THz wavelength range. 

In this paper we present complete iris line theory calculations. In particular, 
the iris line eigenmodes are studied. The analysis of the transmission line has 
been performed in two steps: the first one consists of numerical simulations 
employing the method by Fox and Li [6], the second one by analytical 
calculations using the method by Vainstein [7, 8]. Numerical simulation 
results are in good agreement with analytical results. 

In order to efficiently couple radiation into the transmission line, it is de- 
sirable to match the spatial pattern of the source radiation to the mode of 
the transmission line. To this end, it is advisable to generate radiation from 
the spent electron beam directly in an iris line with the same parameters a 
and b used in transmission line. In this way, the source generates THz radi- 
ation pulses with a transverse mode that automatically matches the mode 
of the transmission line. We developed a theory supporting this choice of 
THz source. As for the microwave waveguide case, one can use a Green's 
function approach to solve the field equations. According to method by 
Vainstein [7], one can set complex boundary conditions for the field on the 
virtual side surface of the iris line, which are called impedance boundary 
conditions. The problem of an open waveguide excitation is thus reduced to 
that of a closed waveguide. Our consideration is quite general, and can be 
applied to edge radiation sources as well as to undulator radiation sources 
in the presence of an iris line. 

We present a complete design for a THz edge radiation source at the LCLS 
baseline. It includes a 15 m-long electron beam vacuum chamber equipped 
with an iris line, and a 100 m-long transmission line with the same pa- 
rameters. The transmission line, which develops through the access maze, 
presents six 90 degrees turns with plane mirrors at 45 degrees as functional 
components. It is possible to match incident and outgoing radiation without 
extra losses in these irregularities. The proposed setup takes almost no cost 
nor time to be implemented at the LCLS baseline. 

The THz edge radiation source is limited concerning the achievable field 
strength. An optimal expansion strategy for the LCLS THz source would 
include a THz undulator source. Such modern, high power THz source is 
exemplified by devices such as the coherent THz undulator source at FLASH 
at DESY [9, 10]. In this paper we will describe this likely extension and its 
accommodation together with the proposed edge radiation source. Both 
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Fig. 1. Optical arrangement for the THz beam transport system, (a) THz beam lens 
system. Focusing of the THz beam can only be performed with reflective optics. 
In fact, lenses made of any material (except diamond) would absorb almost all 
radiation, (b) The use of an iris line, made of periodically spaced metallic screens 
with holes is a way to obtain a high transfer efficiency at small diameter of the 
beam pipe. 

options have advantages, and disadvantages. The edge radiation source is 
characterized by low cost and can operate at all electron beam energies, but is 
limited in field strength. The THz undulator source would provide increase 
in field strength, but a trade-off must be reached between the far-infrared 
and the X-ray achievable frequency range due to technical limitations on 
magnetic field strength of room-temperature electromagnetic undulators 2 . 



2 Principles of THz radiation generation based on the use of an iris line 

2.1 Coherent THz edge radiation source 

As discussed above, the availability of a THz source at XFEL facilities should 
be complemented by the availability of a suitable THz beam transport sys- 
tem, which must guide the beam for distances in the 100 meters range. 

The THz beam can only be transmitted with quasi-optical techniques. In 

2 On the one hand, a decreased electron energy extends the far-infrared frequency 
range. On the other hand, it limits the X-ray frequency range. 
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particular, the idea of providing a periodic phase correction for the free- 
space beam, in order to compensate for its divergence, is very natural. In 
the 1960s numerous attempts of designing various quasi-optical transmis- 
sion lines were reported. In particular, it was proposed to use open beam 
waveguides such as lens guides, mirror guides, and iris guides [6]-[17]. The 
competition among different proposals ended with the victory of mirror 
guides, still in use today e.g. for plasma heating [16, 17]. 

Focusing of the THz beam can only be provided with reflective optics since 
lenses made of any material would reflect and absorb all radiation at long 
distances 3 . Fig. 1 (a) shows the optical arrangement of a transport system 
based on the use of a mirror guide. Each focusing unit is composed by 
two matched copper mirrors. The mirrors are separated in such a way that 
the incident angle is sufficiently small to minimize astigmatism. Existing 
mirror guides are characterized by a maximal length of about 40 - 60 m 
[18, 19]. However, at XFEL facilities, scientists need a beam transport system 
working for significantly longer distances. The transverse size of the line is 
critical for this kind of applications of an open waveguide. Since at LCLS 
baseline the open guide should be installed within a maze with limiting 
size, advantages of the mirror guide are not evident. 

In order to keep the transverse size of the guide to an acceptable level we 
propose, as an alternative to the mirror line solution, to use the iris line in Fig. 
1 (b). Iris lines are characterized by a low attenuation of the fundamental 
mode, self-filtering of higher order modes, wide operational wavelength 
range and mechanical integrity of the structure. The first investigation of 
the influence of diffraction effects on the formation of the field eigenmodes in 
an iris line was carried out by Fox and Li using physical optics techniques 
[6]. According to their explanation, an iris line operates in the following 
way. Consider an electromagnetic wave passing inside the sequence of iris. 
When the wave diffracts at the first iris it produces a diffraction pattern 
in the plane of the next iris. If the second iris lets the main maximum of 
the diffraction pattern through, diffraction losses are minimal. Further on, 
sideband maxima of the diffraction pattern produced by the second iris are 
smaller than those in the first pattern. After the wave has passed a large 
number of irises, a field eigenmode is formed which has low diffraction 
losses. 

A very different and mathematically solid approach to the same problem 
was introduced by Vainstein [7, 8]. His studies were based on direct solution 
of Maxwell equations. They showed that the way the electromagnetic field 



3 In additional to absorption losses, a dielectric lens suffers from reflections at its 
boundaries as well. These surface reflections include a number of effects that are 
almost always harmful. 
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is confined inside an iris line, is essentially different from the guiding mech- 
anism in stable lens guides or mirror guides. In fact, for iris guides, both 
diffraction and reflection from the iris edges are involved. An analysis of 
this effect allowed Vainstein to derive for the first time analytic expressions 
for field distribution and mode losses in iris guides. Complex boundary 
conditions are set on the virtual side surface of the iris line, the so called 
impedance boundary conditions, and the problem of an open waveguide 
transmission is therefore reduced to that of a closed waveguide. An iris line 
has distinct transverse modes. These eigenmodes are similar albeit not iden- 
tical to microwave waveguide modes. For instance, one can find that iris 
line modes, in contrast with microwave waveguide modes, are independent 
of the polarization of the radiation. 

An infinite discrete set of eigenfunctions with corresponding complex eigen- 
values can be found. This set of modes comprises leaky modes which are 
non-orthogonal and non-normalizable with respect to the usual definition 
of inner product. To overcome this difficulty one can define a bilinear form 
which, for our purpose, is equivalent to an inner product under which 
modes are orthonormal. Bi-orthogonality is often exploited in different 
problems [11]. For example, it can be shown that the treatment of a mi- 
crowave overmoded waveguide with resistive walls can be treated in terms 
of Leontovich impedance boundary conditions and of a bi-orthogonal set of 
eigenmodes, in formal analogy with those for a microwave waveguide with 
perfectly conducting walls. From the mathematical viewpoint, the iris guide 
theory with Vainstein boundary conditions is not fundamentally different 
from the theory of microwave waveguide with resistive walls. 

As already mentioned, an iris line constitutes a suitable THz beam transport 
system for the LCLS baseline. Calculations indicate that the losses of the 
principal mode should be in the 10% level per 100 meters. The transfer 
line pipe radius is estimated to be equal to 7 cm. It is therefore technically 
feasible to install such transfer line inside the access maze. However, in order 
to efficiently couple radiation into transmission line, one also needs to match 
the spatial pattern of the source radiation to the mode of transmission line. 
First we analyze a coherent edge radiation source in the THz range, which 
is relatively simple to implement at XFELs. A setup where edge radiation is 
formed with the help of upstream and downstream metallic screens (foils) 
is shown in Fig. 2. The edge radiation from the upstream screen is extracted 
by the downstream metallic screen, which acts as a mirror, and is sent to the 
iris transmission line. The length of the straight section between upstream 
screen and mirror plays the role of the length of the insertion device for 
edge radiation 4 [20]. The matching problem is easily solved if the THz 



4 A hole may be present or not in the edge radiation screens. In the hard X-ray 
regime of operation Be foils without hole can be used. In the soft X-ray regime the 
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in this scheme edge radiation is automatically matched with transport line 



Fig. 2. Scheme for generating edge radiation in an iris line. In this scheme edge 
radiation is automatically matched to the transport line. 

source is equipped with an iris line as well. We consider an electron beam 
moving along the z-axis inside the axisymmetric iris line, Fig. 2. Such source 
generates THz radiation with a fundamental mode that is automatically 
matched to that of the transmission line. We developed a theory of such 
kind of THz source using Vainstein impedance boundary conditions. As 
in the case of a microwave waveguide with resistive walls, one can use a 
Green's function approach to solve the field equations. Only non-azimuthal 
symmetric modes turn out to be driven by the uniform motion of the space 
charge distribution in an axisymmetric iris guide. From a physical viewpoint 
this is sound result. In fact, radiation is related with energy change of the 
particles, which happens through the scalar product of the electric field and 
velocity of the particles. Since the transverse velocity is equal to zero in the 
edge radiation case, symmetric modes cannot lead to any energy change of 
the electrons moving along the axis of the axisymmetric iris guide. Let us 
focus on the fundamental non symmetric mode only. Neglecting losses, for 
the moment, the amplitude for the orthogonal polarization components of 
the field in a cartesian coordinate system, where x and y are the transverse 
horizontal and vertical directions is given by: 



A x = Aicos(0)/i(3.83r/fl) 

A y =A 1 sin(0)/ 1 (3.83r/a), (1) 

minimal size of the hole is defined by the condition that losses of the soft X-ray 
radiation from the baseline undulator due to the aperture limitation should be 
avoided. In any case, a hole with a diameter of a few mm will not perturb the soft 
X-ray beam nor the electron beam, nor the THz beam. 
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Fig. 3. For some experiments a smooth THz field is essential. In these cases, a THz 
bandpass filter should be used. 

where }\ indicates the first order Bessel function of the first kind, a is the iris 
radius, A\ is a constant for a certain longitudinal position, and r = -\Jx 2 + y 2 , 
(p = arctan(y/x) designate polar coordinates. It follows from Eq. (1) that 
the direction of the electric field is radial i.e. varies as a function of the 
transverse position. Therefore, for the edge radiation case we have two 
separate amplitudes for two orthogonal polarization directions. They are 
not azimuthal symmetric because they depend, respectively, on x/r and y/r, 
i.e. on the cosine and on the sine of the azimuthal angle. Only if one sums 
up the intensity patterns referring to the two polarization components one 
obtains the azimuthal symmetric intensity distribution 



I = A\]{ (3.83r/fl) . 



(2) 



Edge radiation is characterized by broad spectrum and radial polarization. 
Many practical applications require the control of the bandwidth and of 
the polarization of the THz pulse. The formation of a THz edge radiation 
pulse usually involves monochromatization and polarization filtering by 
inserting a THz band pass filter for reducing the spectral bandwidth and a 
wire grid polarizer for producing linearly polarized THz radiation, Fig. 3. 

When dealing with a setup where THz radiation from an ultra relativistic 
electron beam is extracted by a mirror (in our case, the Be foil) and sent to 
a THz diagnostic station or user hutch, one usually talks about Backward 
Transition Radiation (BTR). The main problem to solve is in the specification 
of the electric field distribution at some position where the mirror is present. 
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It should be stressed that specification of the field at the mirror position must 
be considered as the first step to the specification of the field at the sample 
position. Such first step is considered separately, because the field at the 
mirror position is independent of the type of mirror and outcoupling optics. 
Once the field at the mirror position is known, the problem of specification 
of the field at the sample position can be solved with the help of Physical 
Optics techniques. 

Let us discuss the problem of field characterization at the mirror position 
in more detail. When electrons are in unbounded space and come from 
an infinitely long straight line, the field distribution in the mirror plane 
can be calculated analytically following Ginzburg and Frank [21]. In most 
practical cases, however, the Ginzburg-Frank equation is not applicable 
because two basic assumptions of the analytical derivation are not fulfilled: 
the electron beam is moving inside the metallic vacuum chamber, which 
acts effectively like an overmoded waveguide, and the straight line has a 
finite length with a bending magnet at the upstream end. It is known (see 
e.g. [20]) that the Ginzburg-Frank theory is the limiting case of the more 
general theory of Edge Radiation (ER) in unbounded space. Emission of 
edge radiation in the presence of metallic boundaries has been a much 
less-treated subject in literature, compared to the unbounded space case. 
To the best of our knowledge, there is only one article reporting on edge 
radiation from electrons in a homogeneous metallic overmoded waveguide, 
in particular with circular cross-section [20]. The method described in [20] 
is therefore capable of treating realistic experimental setup. Here we apply 
method [20] to the case of an iris guide. 

To fix ideas we focus our attention on the setup in Fig. 2. Electrons travel 
through the usual edge radiation setup, similarly as in [20]. The difference 
is that now we account for the presence of an iris guide along the straight 
section. Since electrons pass through an upstream edge screen, one may as- 
sume that the iris guide starts at the upstream screen position. In Section 5 
we will calculate the field distribution at the mirror position, which should 
be subsequently propagated to the experimental hall. The presence of an 
upstream edge screen seems at first glance not necessary, because electrons 
come in any case from the straight section. However, due to variation of 
vacuum chamber cross-section upstream and downstream of the baseline 
undulator, and presence of the undulator magnetic field, characterization 
of the field distribution at the upstream open end of the iris guide is prob- 
lematic. In spite of this, the electric field of the electron beam in the plane 
immediately behind the upstream screen can be well defined as zero, lead- 
ing to extra simplifications. In this case we deal with a well defined problem 
and this allows us to characterize the field distribution at the end of the iris 
guide in the mirror plane. 
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Fig. 4. Spectral properties of undulator radiation. A single electron passing an 
undulator radiates an electromagnetic wave with N w cycles. For on-axis radiation 
the relative spectral FWHM bandwidth is 0.89 /N w near the central frequency. 
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Fig. 5. Some typical spectra of the longitudinal electron distribution and undulator 
spectral distribution function. In general, the bunch form factor varies much more 
slowly in co than the sharp undulator resonance term, and we can replace the form 
factor function by the constant value at the center of the sharp resonance curve. 

2.2 Possible extension plans towards a coherent THz undulator source 

The THz edge radiation source is limited as concerns the maximal field 
strength achievable. The optimal expansion strategy for the LCLS THz 
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source should include a THz undulator source to be installed at the LCLS 
baseline. 

At resonance, undulator radiation can be described by three parameters: 
period A w , undulator parameter K, and number of periods, N w . A single 
electron passing an undulator radiates an electromagnetic wave with N w 
cycles, Fig. 4. For the radiation within the cone of half angle 



the relative spectral FWHM bandwidth is Aco/co = 0.89/Na, near the central 
frequency 



where y is relativistic factor, Fig. 4. The proposed THz undulator is an 
inexpensive, planar electromagnetic device with ten periods, each 1.5 m 
long. At the operation wavelength of the THz source around 0.1 mm, the 
peak value of the magnetic field is about 1 T at an electron beam energy 
around 5 GeV. 



In the case of the THz undulator at the LCLS, the electron beam transverse 
size is much smaller than the diffraction size. This means that, as pertains 
the characterization of the THz pulses, the electron beam can be modeled 
as a filament beam. In this case the electron beam current is made up of 
moving electrons randomly arriving at the entrance of the undulator 



I(t) = (-e)Y j 5(t-t k ), (5) 

k=\ 

where <5(-) is the Dirac delta function, (—e) is the electron charge, N e is number 
of electrons in a bunch, and fjt is the random arrival time of the electrons 
at the undulator entrance. The electron bunch profile is described by the 
profile function F(t). F(t)dt represents the probability of finding an electron 
between time t and time t + At. The beam current averaged over an ensemble 
of bunches can then be written in the form: 



</(*)> = (-e)N e F(t) . (6) 

The radiation power at frequency co, averaged over an ensemble, is given 
by the expression: 
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(P(co)) = p(co)[N e + N e (N e - l)\F(cv)\ 2 ] , 



(7) 



where p(co) is the radiation power from one electron and F(co) is the Fourier 
transform of the bunch profile function 5 . For wavelengths shorter than 
the bunch length the form factor reduces to zero. For wavelengths longer 
than the bunch length it approaches unity. A sample of undulator radiation 
spectrum is shown in Fig. 5. The distribution of the radiation energy within 
different harmonics depends on the value of the undulator parameter K. In 
our case of interest K » 1. Strong undulator maxima are present at co = 3co , 
co = 5co and so forth. The energy measured by the detector is proportional 
to the convolution of the square modulus of the bunch form factor and of 
the spectral line of the undulator spectrum, as is illustrated in Fig. 5. Note 
that the bunch form factor behaves like an exponential function at high 
frequencies, and falls off rapidly for wavelengths shorter than the effective 
bunch length. For wavelength about three times shorter than the effective 
bunch length, the radiation power is reduced to about one percent of the 
maximum pulse energy at the fundamental harmonic. As a consequence, 
sharp changes of the bunch form factor result in attenuation of the higher 
undulator harmonics. In general, |F(cc)| 2 varies much slower in frequency 
than the sharp resonance term, and we can replace |F(cu)| 2 by the constant 
value |F(cuo)| 2 at the center of the sharp resonance curve. 

Let us consider the practically important case of an electron bunch with 
strongly non-Gaussian shape similar to that used to drive the LCLS. Fig. 6 
shows the current distribution along the bunch [22]. The nominal charge is 
0.25 nC. The electron bunch has a complicated shape, which is reflected in 
the squared modulus of the form factor shown in Fig. 7. In the case of the 
LCLS, the squared of the bunch form factor modulus falls off rapidly for 
wavelengths shorter than 0.06 mm. At the opposite extreme, the dependence 
of the form factor on the exact shape of the electron bunch is rather weak 
and can be ignored for wavelengths longer than 0.1 mm. 

Similarly as for the edge radiation case, in order to solve the matching 
problem with the transmission line, we propose to generate the coherent 
THz undulator radiation directly in the iris line, Fig. 8. 

Summing up, in view of the practical application to the THz undulator line 
at LCLS, there is a need to develop a comprehensive theory of undulator 
radiation in the presence of an iris line. As for the edge radiation source 
case, one can use a Green's function approach to solve the field equations, 
and the paraxial (overmoded waveguide) approximation can be used. In 

5 A similar expression holds for the case of edge radiation. In fact, the electron 
beam transverse size is much smaller then the THz edge radiation diffraction size. 
This means that, as pertains the characterization of the THz edge radiation pulse, 
the electron beam can be modeled as a filament beam. 
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Fig. 6. Electron beam current profile at the LCLS, after [22]. 

addition to the paraxial approximation, in the undulator radiation case, the 
resonance approximation can be exploited too. We thus consider a large 
number of undulator periods and a frequency range of interest close to the 
fundamental harmonic. In free-space and under resonance approximation 
the radiation from a planar undulator is horizontally polarized, and thus 
constitutes a replica of the undulator polarization properties. Moreover, the 
field exhibits azimuthal symmetry. These properties are unvaried when an 
axisymmetric iris line is introduced. The wiggling amplitude of the electron 
in the undulator is taken to be small with respect to the dimension of the 
waveguide. This greatly simplifies analytical calculations, and describes our 
practical case of interest. 

Applying the Green's function approach in this case one obtains a space 
invariant polarization and an azimuthal symmetric field distribution. In 
an axisymmetric iris guide only azimuthal symmetric modes turn out to 
be driven by the wiggling electron moving along the z axis. Let us focus 
on the first, dominant azimuthal symmetric mode. Neglecting losses, for 
the undulator radiation case we expect the following amplitude for the 
horizontally polarized radiation pulse: 

A x =A J Q (2Ar/a), (8) 
where Jo is the Bessel function of the first kind of zero order, Fig. 8. 
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Fig. 7. Squared modulus of the electron beam form factor, \F\ 2 , corresponding to the 
current profile in Fig. 6. 
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in this scheme THz undulator radiation is automatically matched with transport line 



Fig. 8. Scheme for generating THz undulator radiation in an iris line. The THz 
undulator radiation is automatically matched to the iris transport line. 

In contrast to edge radiation sources, THz undulator sources operate within 
the typical spectral window of undulator radiation Aco/co ~ 1/N W and the 
radiation is almost completely linearly polarized, Fig. 4. Additionally we 
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Fig. 9. A practical arrangement for the coherent THz source focusing system. 



should consider the physics of the focusing process. In fact, at the end of 
the transmission line the THz pulses should be focused with the help of 
a converging lens or a paraboloidal focusing mirror, Fig. 9. The focusing 
spot of the first symmetric mode of an iris guide should be compared with 
the focusing spot of the first non symmetric mode of the same iris guide. 
Summing up, in view of practical applications, an important advantage of 
THz undulator sources over edge radiation sources is a few ten times (for 
N w ~ 10) higher intensity on the sample as compared to THz edge radiation 
sources at the same electron beam parameters. 



3 Theoretical background of an iris line 



We begin our theoretical investigation from Fig. 10, a sketch of an iris 
line, which is formed by perfectly absorbing screens of infinite extent with 
collinear circular holes. The radius of the holes is indicated with a, while the 
distance between two screens is indicated with b. 

An analysis of the behavior of an electromagnetic wave passing through the 
line was given first in terms of numerical calculations by Fox and Li [6], and 
then in terms of analytical results by Vainstein [7, 8]. 
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Fig. 10. Iris line comprising a series of collinear identical apertures cut into parallel 
and equally spaced black (i.e. perfectly absorbing) partitions of infinite extent. 

3 . 1 Numerical method 



Consider the iris guide in Fig. 10, and call E(r', 0) the transverse profile of a 
given polarization component of the slowly varying envelope of the electric 
field on the iris at z = in the space-frequency domain, that is at a given 
fixed frequency co. The field at the iris at z = b is easily obtained as 



E(r,b) = dr'K{r,r';b)E{r',0) , 



(9) 



where function K(r, f; z) is a Green's function, or propagator, for the paraxial 
wave equation which describes the field evolution. In physical optics one 
uses the Fresnel propagator which is the solution of the wave equation for 
unbounded space (with zero boundary at infinity, outgoing wave). In the 
case of a circular shape of the iris one may use polar coordinates and obtain 
[6]: 



2n a 



E(r,(p,b) = j d(p' j dr'r'K(r,r'(p,(P';b)E(r',(p',0) , 



(10) 



o o 



with 
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K(r,r' ( p, ( p';b) = 



ico exp[-ikb] 



cb 



exp —ik 



( 



r 2 + r 



2b 



./2 



yCOS^-^')) , (11) 



which is valid at large Fresnel numbers a 2 /(bA) » 1. Iterations can be per- 
formed iris after iris, to obtain the field at the position of any iris at arbitrary 
distance from the entrance. The modes of the waveguide are, by defini- 
tion, field distributions that do not depend on the longitudinal z coordinate. 
Since we are dealing with the modes of an empty waveguide with no gain 
medium, the mode amplitudes decreases along the z axis due to diffraction 
at the irises. However, the shape or spatial pattern of the field distribution 
defining a given mode does not. It follows that the field profile of a certain 
mode at the n + 1 iris, E(x, y, (n + l)b) is simply E(x, y, nb) multiplied by some 
complex number 6 y. It follows that the following integral equation holds 



where for simplicity we drop the explicit reference to the z dependence. 
Note that Eq. (12) may be written in the operator form 



where T is the operator corresponding to a trip between two successive 
irises: 



According to Eq. (14), the modes of the iris guide are the eigenfunctions 
of the operator Y for the guide. The number y is, instead, the eigenvalue 
corresponding to the eigenf unction E(f). In contrast with usual quantum- 
mechanical situations, where eigenvalue equations include Hermitian op- 
erator, here T is generally not Hermitian. For instance, its eigenvalues are 
complex, whereas the eigenvalues of a Hermitian operator are always real. 

Due to diffractive losses at the irises, the total energy associated with the 
field inside the guide diminishes, and one must have \y\ < 1. The iris guide 
modes may be found by solving Eq. (9) numerically. This is usually done 
according to the method developed by Fox and Li in 1961. One starts by as- 
suming some initial field E(x, y, 0) on an iris, usually just E(x, y, 0) = constant. 
The field is then propagated to the other iris by calculating the integral in 

6 Not to be confused with the relativistic factor y of the electrons. 




(12) 



TE = yE, 



(13) 




(14) 
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Eq. (9) numerically and iterating until the field on the irises is unchanged 
(within some prescribed numerical error) on successive iterations, except 
for a constant factor y. The field yielded by this method is a solution of 
Eq. (12), that is, it is a mode of iris guide. Only one mode, that is the one 
with smaller losses per iteration will be returned following this procedure. 
However, higher-loss modes can also be obtained. For example, looking for 
solutions with the form 



E(r,tp,z) = R„(r,z)exp(-in(p) , 



(15) 



with n = 0,1,2..., Eq. (10) and Eq. (11) can be rewritten for R„(r,z) after 
integrating in d(p' as 



a 

R„(r,b) = J dr'r'K r (r,r';b)R n (r',0), 



(16) 



with 



i n+ co T / krr' 



K r (r,r';b) = — U — lexp 



-tie 



., r 2 + r' 2 



2b 



(17) 



which allows to find different modes as n varies. 



3.2 Analytical method 



The analytical method by Vainstein is based on the introduction of partic- 
ular, complex boundary conditions, which are called impedance boundary 
conditions. Once these boundary conditions for the field are formulated on 
the virtual side surface of the iris line, the original problem, which is about 
an open waveguide excitation, is simplified to that of a closed waveguide. 

In order to understand how the impedance boundary conditions come into 
play, we consider, following [7, 8], the simpler case of diffraction of a plane 
wave at a stack of semi-infinite screens, as depicted in Fig. 11. We assume 
here that the screens are totally absorbing. When the wave illuminates one 
screen, the total field can be presented as the sum of the unperturbed wave 
and a cylindrical wave generated by an image source at the edge of the 
screen. The unperturbed wave propagates up to the following screen and 
produces another cylindrical wave at the second edge. As a result, the 
diffracted wave is given by the sum of many cylindrical waves generated 
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Fig. 11. Diffraction of plane wave at semi-infinite screens. 




Fig. 12. Geometry of the diffracted wave. 

at the edges of the screens. Constructive interference takes place only along 
a series of privileged directions. This is shown in Fig. 12. The reader can 
easily recognize that the path difference between A\BiC\ and A 2 B 2 C 2 is just 
given by D^B^ - D 2 B 2 - b(0 2 - a 2 )/2. Constructive interference is given for 



kb(6 2 - a 2 ) 12 ^ Inn (18) 
with n = 0,1,2..., while it is zero at all other directions. For n = one 
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obtains the zeroth order diffracted wave, which is just "reflected". When 
the angle of incidence is small a 2 «: A/(2nb), the screen does not absorb, but 
rather reflect the wave. The finding of the reflection coefficient pertaining 
the zeroth order diffraction is due to Vainstein, and reads 



Ro = - exp [-j8 (l - i)a(kbf 2 ] , (19) 

with 7 j6 = 0.824. Note that Rq is complex, meaning that the zeroth order 
diffracted wave suffers phase shift and attenuation. The impedance bound- 
ary condition are a consequence of Eq. (19). 

In fact, the slowly varying envelope of the total field in the plane of the 
screen turns out to be written as 



E(x) = A[exp(ik x x) + R exp(-ik x x)] , (20) 

where A is just the amplitude of the incident field. Then, using a k x /k z 
and k 2 +k 2 = 4n 2 /A, and taking advantage of Eq. (19) one obtains that the 
logarithmic derivative of E at the edge of the screen can be written as 



d]n(E) IdE 1-Ro 

— — \ x=0 = =-Hx=o = ik x . (21) 

ax £ ax 1 + Rq 

Assuming a ^Tkb <sc 1 and using the expansion exp[-j6 (l - i)oc Vfcfc] 1 - 
f$o(l - i)a ^Tkb, we obtain 



^T Uo = f ^ lx=0 = " [^ o(1 + V^W^)]" 1 • ( 22 ) 

The main feature of this formula is that the term k x is completely excluded 
from it, which allows one to use Eq. (22) as an approximate boundary 
condition. Note that this is valid not only for plane waves. In fact, any 
wave can be decomposed in terms of a linear superposition of plane waves. 
Since Eq. (22) is valid for each component, it must be valid for their linear 
superposition as well. Eq. (22) can also be extended for any shape of the 
(virtual) boundary surface, provided that the typical value of the curvature 
radius is much larger than the wavelength. Eq. (22) is named after Vainstein, 
who first derived it, and can be written in vector form as 



jSo turns out to be related with one of the most famous mathematical functions, 
the Riemann zeta function C- In fact, it is given by -£(1/2)/ yfn, see [7, 8]. 
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E + (1 + i)j8o y/cb/(4co) (it ■ V ± )E 



= 0, 



(23) 



where S is the cylindrical virtual surface of the iris guide, and ft is the unit 
vector normal to S. Let us consider axially symmetric iris line. We seek a 
solution for the field amplitude E(r, (p, z) in the form 



E = u n j(r) exp[-in<ft - ik z z] , 



(24) 



with n = 0,1, 2, ... . The functions u nk (r) are subjected to the following homo- 
geneous equations 



r 2 < ; . + ru' nj + [(k nj ) 2 - n 2 )]u nj = , 



(25) 



and satisfy the boundary conditions 



[u nj + (1 + 0j8 y/cb/(4co)u' nj ] r=a = . 



(26) 



In the first order of the small parameter M = (8nN) 1 ^ 2 , where N = a 2 /Ab is 
the Fresnel number with a the iris radius, the functions u n j assume the form 



Unj = Jn(k nj r) (27) 
where 



V n] 



knj = —[l-(l+i)PoM\, (28) 



a 

and v n j is the ;-th root of the n-th order Bessel function of the first kind (i.e. 
Jn(Vnj) = 0). Substituting the expression for k n j into the dispersion relation 

^/ = ^ (29) 
we obtain the following expression for k z : 

k z b = ^- 2v 2 nj M 2 + Av 2 nj M 3 (l + OjS . (30) 

For the eigenmode with transverse wavenumber k n j, the fraction of the 
radiation power losses per transit of one iris is given by 
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Fig. 13. Comparison of the fundamental non-symmetric mode derived analytically 
and numerically for the Fresnel number N = a 2 /(Ab) — 20. 



21m(k z b) = 8v 2 nj M 3 (i . 



(31) 



The relative loss of the /-th mode of order n after traveling for a distance z 
is therefore given by 



'AW) =1 

/ r, 



w 



exp 



'a) 



{2nNf 2 b 



= 1 - exp 



(27l) 3 / 2 fl 3 & 



(32) 



Due to the exponential dependence on v 2 ., only the lower order modes 
tend to survive. Note that the exponent in Eq. (32) depends on the distance 
between two irises, b, only weakly as yb, while there is a much stronger 
dependence on A and a. 



3.3 Comparison between analytical and numerical methods 



Comparisons between analytical and numerical results are given in Fig. 13 
and Fig. 14, respectively for the edge radiation and undulator cases. For the 
sake of completeness, Fig. 15 shows a comparison between the phase of the 
main symmetric mode calculated analytically and numerically as well. 
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Fig. 14. Comparison of the fundamental symmetric mode derived analytically and 
numerically for the Fresnel number N = a 2 /(Ab) = 20. 

Note that there are some minor differences between analytical and numerical 
results. These differences, in form of ripples in the numerically calculated 
modes, are not a residual from higher order modes. In fact, they do not 
disappear as the number of iterations increases. They can be explained 
noting that Eq. (22) only accounts for the zeroth order diffraction, see Eq. 
(19) with n = 0, while the numerical method by Fox and Li accounts for all 
diffraction orders. As a result, some difference should be expected between 
the iris line modes calculated with the two methods. 

A comparison between Eq. (32) and losses from numerical calculations is 
given in Fig. 16 as a function of the Fresnel number for n = 0, and in Fig. 
17 for n = 1. The convergence of these results as a function of the iteration 
number is exemplified in Fig. 18 in the case N = 20. After a large number of 
iteration, the power losses between two successive irises becomes constant. 
This constant is the result finally plotted in Fig. 16 and Fig. 17. 

3.4 Technical requirements 

It is important to study the influence of imperfections in the construction of 
the iris line on its nominal properties. With little modification, the computer 
technique by Fox and Li can be used to study the effects of misalignment of 
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Numerical results 

Analysis 



N=20 




Fig. 15. Comparison of the phase of the fundamental symmetric mode derived 
analytically and numerically for the Fresnel number N — a 2 /(Ab) = 20. 

iris positions in the longitudinal direction. The importance of this effect is 
shown in Fig. 19 for A = 0.1 mm, b = 30 cm, and a = 1.35 cm. When the irises 
are adjusted with an accuracy of 1 mm (rms) the irregularities of the iris line 
do not result (with graphical accuracy) in additional diffraction losses. The 
result that the iris line is a "non-resonant" device should not be surprising, 
and can be deduced from the explanations presented above in this Section. 

Let us now study the requirements on the accuracy of the transverse align- 
ment between the irises. When one of the iris is shifted off the axis, it may 
cause distortion in amplitude and phase of the diffracted wave. To estimate 
the change in amplitude, one should remember that the amplitude of the 
wave produced by the edge of the iris (see Fig. 11) in the region of the next 
iris is proportional to the Fresnel integral 



F(i/0 = J exp(zx 2 )dT , (33) 

where \p = x ^n/(Ab). The Fresnel integral F(ip) is decreased by about a 
factor 2 when x changes from to yfXbfn. Therefore, one can consider the 
region for x > ^jAb/n as the region of the shadow. If an iris is shifted off 
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Fig. 16. Power losses per transit between two successive irises vs N = a 2 /(Ab) for 
the first non-symmetric mode. The curve is calculated with the analytical formula 
Eq. (32). Circles are the result of calculations with the Fox-Li iterative method. 



the axis by a value ~ ^Ab/n, either such iris or the next one fall into the 
region of the shadow and does not produce any diffracted wave. To estimate 
phase errors, we consider Fig. 12. One can see that a shift of the iris by a 
certain value bx in the transverse direction causes a phase shift of the mirror- 
reflected wave of about laco ■ 5x/c. A phase shift equal to n is achieved at 
5x = A/ (4a). In paraxial approximation a 2 <sc A/(nb), and the value of the 
admissible shift of the iris is mainly defined by aperture restrictions, and 
not by phase distortions. Summing up, requirements on the accuracy of the 
transverse alignment of the irises is given by 5x <^c s/Ab/n. For A - 0.1 mm, 
b — 30 cm, we estimate a restriction 5x «: 3 mm. 

Finally, with little modifications, the computer technique employed up to 
now can be used to study the effects of particular aperture restrictions. 
When the radius of the irises are manufactured with an accuracy 5a, these 
irregularities result in extra diffraction losses. The importance of this effect is 
shown in Fig. 20. Simulations show that when the irises are manufactured 
with an accuracy better than 1 mm in both transverse and longitudinal 
direction, manufacturing errors do not result in extra diffraction losses in 
our case of interest. 



In conclusion it should be noticed that requirements on the transverse extent 
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Fig. 17. Power losses per transit between two successive irises vs N = a 2 /(Ab) for 
the first symmetric mode. The curve is calculated with the analytical formula Eq. 
(32). Circles are the result of calculations with the Fox-Li iterative method. 
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Fig. 18. Power losses between two successive irises at N - 20 as a function of the 
iteration number. Left plot: first symmetric mode. Right plot: first antisymmetric 
mode. 



of the screens i.e on the difference between the pipe radius R and the hole 
radius a is given by R - a » ^Ab/n 3 mm. Additionally, there is no need to 
use completely absorbing screens in the iris line. It would even be preferable 
to use reflecting screens. In fact, within the accepted limitations, reflecting 
screens are almost identical to absorbing screens as concerns diffraction 
effects. 
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Fig. 19. Power of the principal symmetric mode versus number of transits. Cal- 
culations have been performed with the Fox and Li method. Here the wavelength 
A = 0.1 mm, the iris separation is b — 30 cm, the hole radius is a - 1.35 cm. The solid 
curve corresponds to the perfect iris line. The circles represent results from numer- 
ical simulations when distances between irises are randomly distributed around 
b — 30 cm with rms o\, - 1 mm. The inset is an enlargement of a region of the main 
plot. 




Number of transits Number of transits 

Fig. 20. Simulated power losses due to aperture distortions as a function of the 
number of transits. The radia of irises are randomly distributed with an rms value 
of 1 mm (left plot) and 0.5 mm (right plot). Here the average hole radius is a — 1.35 
cm, b =30 cm, A = 0.1 mm. 

3.5 Focusing of fundamental edge radiation and undulator radiation modes 

For practical experimental purposes, edge or undulator THz radiation needs 
to be focused onto the sample. This can be done by placing a focusing mirror 
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Fig. 21. Single lens focusing system for the edge radiation source. 
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Fig. 22. Single lens focusing system for the undulator radiation source. 

at the exit of the iris guide, and by positioning the sample in the focus. We 
will treat the mirror as a thin lens with focal length /. The focusing system 
for the edge radiation and for the undulator radiation sources are illustrated 
in Fig. 21 and Fig. 22 respectively. 

In general, in the space-frequency domain, the relation between the slowly 

-> 

varying envelope of the field E(0, f) at the position immediately in front of 

-> 

the lens, Zi ens = 0, and the slowly varying envelope of the field E(f, rf) at the 
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sample position z sam pie = // is given by 



ico 



E(/ ' r ^ ) = 2^7 exp 



icon 



2cf 



! 



icovf ■ r' 



(34) 



where P(r') indicates the pupil function, and will be assumed here to be 



P(r) = 1 for r < a, P(r) = otherwise . 



(35) 



Differences between the edge radiation and the undulator radiation case 

-> 

arise due to the different expressions for E(0, r'). In particular, the radial 
polarization of the edge radiation pulse leads to a non-azimuthal symmetric 
response, whereas the linear polarization of the undulator radiation pulse 
yields an azimuthal symmetric response. 



From Eq. (1) we know that in the edge radiation case 



E(0, f) = A x h p^-\ [cos($'% + sm{cb%] > 



(36) 



while from Eq. (8) we know that in the undulator case the field at the position 
immediately in front of the lens is 



E(0,r')=A / c 



)e x . 

a I 



(37) 



Substitution of, respectively, Eq. (36) and Eq. (37) in Eq. (34) yields the 
following results for the field at the sample in the case of an edge radiation 
pulse 



E(f,?f) = 



coA\ 
cfa 2 



li 



rjco 2 a 2 ^ 1 



c 2 f 



vii/o(vn) exp 



icon 



2cf 



coa \ 

h W>)' 



(38) 



and of an undulator pulse 



icoA c 



2 2 2\ - l 

rjco z a z N 

2 fl 



c 2 f 



Voi/i(v i)exp 



icon 



2cf 



coa \ 



(39) 



It is possible to give a quantitative measure of how well one can focus edge 
radiation pulses compared to undulator pulses. This measure is simply 
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given by the ratio between the maximum of intensity in the two cases, 
assuming that the energy per mode per unit spectral interval is the same. 

First note that the energy per unit spectral interval per unit surface after 
focusing is proportional to the squared modulus of Eq. (38) and Eq. (39), 
respectively for the edge radiation and for the undulator case. The assump- 
tion of same energy per mode per unit spectral interval can be enforced by 
means of the following relation between the field amplitudes: 



(dlf = (40) 

where A and A 1 enter in the expression for the field as in Eq. (39) and Eq. 
(38). The energy per unit spectral interval per unit surface is given by 



dW c 
dcodS 47i 2 



(41) 



Note that the transverse dependence of Eq. (38) and Eq. (39) is in the nor- 
malized distance ft = TfCoa/(cf). For the non-symmetric mode case, the 
maximum of energy density in the focal plane occurs at ff ^ 2.84. For the 
symmetric mode case, such maximum occurs at ff = 0. The relative effi- 
ciency with which the non-symmetric mode can be focused compared with 
the symmetric mode is given by r\ ~ 1. The outcome of this calculation is 
different from what one would have intuitively guessed. If we assume the 
same energy in symmetric and non-symmetric modes, after focusing we 
obtain practically the same maximum energy density on the sample. Our 
consideration is only valid for the aberration-free lens case and within a 
paraxial treatment i.e. for numerical aperture of the lens NA = a/f < 0.3. 
In fact, the field distribution of the non-symmetric mode, Eq. (36) is quite 
different compared to the symmetric mode distribution Eq. (37), and the 
aberrations can have a significant effect on the relative mode efficiency r\. 



4 Green's function for an axisymmetric iris line 



Let us consider the motion of an ultrarelativistic electron in an axisymmetric 

iris line, and indicate with E(r,z,co) the Fourier transform of the transverse 
electric field generated by the electron. Since the radiation formation length 
is much longer than the reduced wavelength K = c/co, the electric field en- 

velope E = E exp [-icoz/c] turns out to be a slowly varying function of z with 
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respect to the wavelength, and it obeys the paraxial Maxwell equation 1 



DE(z, r) = g(z, r) 
The differential operator D in Eq. (42) is defined by 



(42) 



2 2 /a i d 



c dz t 



(43) 



where V ± 2 is the Laplacian operator over transverse cartesian coordinates. 
The vector g(z, r) is specified by the trajectory of the source electron, r'(z), 
and is written as 



8 = 



4ne 



exp 



i ^ dz- 



to 



2y\(z)c 



—v(z) - V ± 

c z 



6 (f- r'(zfj . 



(44) 



Here we substituted v z (z) with c, based on the fact that 1/yl « 1. Eq. (42) is 
Maxwell's equation in paraxial approximation. 

The boundary conditions are taken as in Vainstein [7, 8] as Eq. (23). Since Eq. 
(23) is valid separately for each component, we can solve Maxwell equations 
separately for the field along the horizontal and the vertical directions. Inde- 
pendently of the actual shape of g, the full problem for a certain polarization 
component can be written as 



/ 



DE(z, r) = g{z, r) 

[E + (1 + 0j8 y[cmw) (ft ■ V ± )E] s = . 



(45) 



We proceed by using a Laplace transform technique, that allows one to 
dispose of the partial derivative with respect to z in D in the first equation 
in (45). 

First we define Laplace X and inverse Laplace X -1 transform of a function 

f(z, r) as 

oo 

/(p, r) = £ [/(z, f)] (p) = J f(z, r) exp [-pz] dz (46) 

o 

8 Here and everywhere in this paper we make consistent use of Gaussian units. 
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with Re[p] > 0, and 



6+ioo 



f(z,r) = £- 1 [f(p,i)](z)= J f(p,i)exp\pz]dp, (47) 



6-ioo 



where 5 is a real number larger than all the real parts of the singularities of 

m 

Consistently applying a Laplace transformation to the equation set (45) we 
obtain a redefinition of the problem in terms of the Laplace transform of the 

field E(p,i)=£\E(z,T)](p): 

DE(p, r) = g(p, t) + 2 -fE(0, t) 
[£ + (l+Oj3oV^M(n-Vx)£] s = ' 



where g is the Laplace transform of g and 



2) = (V ± 2 + ^]. (49) 



Note that the presence of the initial condition E(0, r) in Eq. (48) refers to the 
possibility of introducing an external field into the system. In what follows 9 
we set E(0, r) = 0. 

Now, suppose that we find a scalar function G such that 



= J G{r,r',p)g{f,p) dr' . 



(50) 



In this case, the inverse Laplace transform of G, that will be simply written 
as G, is the Green's function for the problem (45), inclusive of the proper 
boundary conditions. As a result 



z 

E(r, z) = j dz' j df G (r, r', z-z')g {f, z') , (51) 



o 



9 In our case of interest we have an edge radiation screen at at z = 0, and at the 
position immediately behind the screen the field of an electron is equal to zero. 
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where we integrate up to z' = z because the radiation formation length for 
z-z' < is very short with respect to the case z-z' > 0. Summing up, we 
first have to find G, then to apply a Laplace inverse transform in order to 
get G and, finally, to solve for E. 

We start by specifying the eigenvalue problem associated with the problem 
set (48), that is 

;w 11 (52) 

[Fj + (1 + Oft V^7(4^) (it ■ V ± )Fj\ s = " 

Setting by definition Ay = 2icop/c- Aj yields the following, equivalent eigen- 
value problem: 

\lF j (r) + A j F j ( i) = 
[Fj + (1 + Oft y[cmw) (it • V ± )F ; ] s = ' 



Boundary conditions are homogeneous, so that the domain of the Laplacian 
operator is the vector space of twice differentiable (square integrable) func- 
tions obeying boundary conditions in (53). However, the Laplacian operator 
defined in this way is not self-adjoint with respect to the usual inner product 
definition, < f,g >= f s drf*g. This is a result of the fact that the boundary 
condition in Eq. (53) are not real. Then, eigenvalues are not real, nor eigen- 
functions are orthogonal with respect to the usual inner product definition. 
In general we do not know wether the spectrum is discrete, completeness 
is not granted and we cannot prove the existence of a set of eigenfunctions 
either. Yet, direct calculations show that 



(FjrF^ J F r Fid^=b fi . 



(54) 



It follows that functions Fj form a bi-orthogonal set of eigenfunctions. This 
allows us to decompose G as 

- F i F i 

G = Y W- 7 ± ^ > ( 55 ) 

2icop c -Aj 



similarly as for a Sturm-Liouville problem. Note that now eigenvalues and 
eigenfunctions are complex. Also note that the dependence of G on p is at 
the denominator in Eq. (55). As a result, G goes to zero uniformly as p — > oo. 
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It follows that the Mellin transform of G, that is the Green's function G, can 
be found with the help of Jordan's lemma and reads: 



G = £jRes[Gexp(p ; z)]} 



(56) 



where 



(Res [G exp[ P; z]]) = ^ Jim ^ |(p + ^) G(p) exp[p^ 



2z'cu 



F ; (r)F 7 (r")exp 



2^7 



One therefore obtains the following expression for G: 



(57) 



' 2o7 



(58) 



Given the symmetry of the problem we now introduce polar coordinates 
(r, <p), define A = (1 + /')/3 ^cb/^coa 1 ) and rewrite Eq. (53) as 



(& + ii + b&)W,<t>) + A j F,ir,<l>) = 
[F 7+ flA<9 r F ; L = 



(59) 



Posing ; = {n, k}, we look for solutions of Eq. (59) of the form 



Fj^nM = f nk exp(-in(p) , 
yielding naturally 



(60) 



F MnM = A nk f n ( VX^r) exp(-in(p) , 



with 



(61) 



(62) 



where }' n indicates the first derivative of the Bessel function of the first kind 
of order, /„, with respect to the argument. Assuming a large Fresnel number 
N = a 1 1 (Kb) » 1, which also implies |A| «: 1, one can solve Eq. (62) for A nk 
and find 
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A nk = -f(l - A) 2 - -f (1 - 2A) , (63) 

where v nk is the k-th zero of /„, that is J„(v nk ) = for fc = 1, 2, 3... . 

Normalization according to the product in Eq. (54), < F ir Fj >= 6 !; finally 
yields 



™ 2 / 2 +1 (v w) 

The following, final expression for the Green's function G follows: 



fc(l-2A) yi yi exp[-m((/) -(/)')] / ic(z - zQv^(l - 2A) \ 
= " W itfe /Li(^) 2ca 2 j 

x/„(^(l-A)r)/„(^(l-A)r'). (65) 

It can be shown that Eq. (65) is in agreement with the findings in [23], where 
a decomposition of the field in azimuthal components is proposed, and 
a Green's function is obtained for each azimuthal component following a 
different reasoning. Eq. (65) can now be used to find the slowly varying 
amplitude of the field along any given polarization component according 
toEq. (51). 

In order to verify the correctness of Eq. (65) we study the free-space limit, 
that corresponds to the limit for large values of the iris radius a — > oo. 

Since we are interested in characterizing the fields over a finite transverse 
direction, and since sources have a finite transverse size, the limit a — > 
oo allows one to substitute Bessel functions in Eq. (65) with asymptotic 
expressions for k » 1. First remember that 



/„(Q* a/4cos(c-^-t), - ' 1 • (66) 



One sees that 



7i£ V 2 4 



Vnk = n(^+k+-^ , (67) 
yielding 
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J 2 n + l( V «k) ~ 



nv nk 



Substituting Eq. (68) in Eq. (65), posing E, = nk/a and replacing the sur 
k in Eq. (65) with an integral over dE we obtain: 



oo 

oo p 

" = -°° o 

x/„(^)/ n (^')exp[-m((^ -<£')] . 
Eq. (69) can also be written as 



ic(z-z') 2 
2co 



oo 

G (r,r',(p-(p\z-z') = J dEEexp 



2ZCU7T 



oo 

oo „ 



dEE exp 



ic(z-z') 2 
2a; 



zc(z-z') 2 
2w 4 



x/ n (^)/ n (^)cos[n((/) -<£')] • 
Here z > z'. The integrals in dE can be performed yielding 



G (r, r , (p - cp ,z - z ) = - A ^_^ Jo f 



47i(z-z') \c(z-z') 



x exp 
x exp 



z'(r 2 + r /2 )cu 
2c(z - z') 
i(r 2 + r' 2 )co 



corf 



2c(z - z') 
Using the Anger-Jacobi expansion 



- — Yr n T , 

2n*-L Jn \c(z-z') l 

n=l ^ v 'I 

cos[n({/> - cp')] . 



HQ + 2 £ r n /„(Q cos[n(<f> - <£')] = expHC cos(4> - <£')] 



we find that in free-space the Green's function Eq. (65) reduces to 



G = —- — : r exp 

4tz(z - z') r 

as it must be. 



zcu 



I -* y r y 

2c(z - z') 
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5 Edge radiation in an iris line 



Let us consider the case of an electron emitting edge radiation in an iris 
line. The vector g, whose horizontal and vertical components are needed in 
order to find the field according to Eq. (51) can be written as 

*=-— exp M ± (74) 

After integrating Eq. (51) by parts, the horizontal and vertical polarization 

-> 

components of E turn out to be given by 



L/2 

E x ,y = ^ J dz'exp^) fdt>d* y G(?,i>,z-zr)6(i>), (75) 

-L/2 

where the choice of the integration limits implies that the reference system 
now has its origin in the center of the drift between the two edges. Since in 
Eq. (65) the Green's function is expressed in terms of polar coordinates, we 
express likewise the partial derivatives with respect to x' and y' as 



sin ((/)') 

d x > = cos((p)d r > — dp (76) 



and 



cos(d)') 

d y , = sm(<t>')d r ^—dp . (77) 

Note that one can transform the partial derivative with respect to y' to that 
with respect to x' by just substituting cos{(\)') with sin((/)'), and sin(^') with 
- cos((/)'). Remembering that 5(f) = 5(r')/(nr'), using 
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exp(-in(p') cos((p') = n for n = ±1 , and otherwise , (78) 

o 
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^ d(p' exp(-in(p') sin((p') = ±n/i for n = ±1 , and otherwise , (79) 

o 
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^ fa t (l-AK\ lim ^(^)_,v, t (l-A) 

11:1 I fl / 



= ±- 



2a 



(80) 



and remembering that v nk = v_ nk , ]^{v\i) = / 2 (v_ 1)t ) and ]\(x) = one 
finds the following expression for the field at the mirror position: 



E x (r,(p) = 



2ie(l - A)(l - 2A)Lcos(<t>) v lfc Mv lfc (l - A)r/a) 



CO 



L 

k=\ 



xsmc 



co_ cv 2 ,(l-2A) W 
2 + 2coa 2 2 



x 2cy 



expf 



4( 



fcv 2 ,(l-2A) L ^ 



2o«? 2 



for the horizontal component and 



Ey(r,<t>) = 



2ie(l - A)(l - 2A)Lsin(<p) ^ v lk h(v lk (l - A)r/a) 



CO 



E 



xsmc 



CO 



2cy 2 



+ 



cv 2 fc (l-2A) W 
2coa 2 2 



k=l 



exp| 



icv 2 (I -2A) L \ 



2coa 2 



for the vertical one. 



(81) 



(82) 



Eq. (81) and Eq. (82) can be used as a basis to calculate the energy per 
unit spectral interval per unit surface at the sample position. In our case 
of interest, the distance between sample and extracting mirror, L s , is much 
longer compared with the length of the edge radiation setup. Assuming 
L s » L and accounting for coherence through the form factor F(co) the field 
energy per unit spectral interval per unit surface at the sample position is 
given by 



dW 
dcodS 



4ti 2 
4e 2 L 2 c 



N 2 e \P(co)\ 
N 2 



4n 2 a 6 co 2 ^ e 
L 



k=i h 



Ik 



xsmc 



co v 2 lk c\ 
4 \cy 2 coa 2 j 



(Vlfc) 

*(?)■ 



ex P -TJT^ L s I 



coV 2 a 3 



(83) 



It should be noted that in the calculation of the square modulus of the 

field, E , crossed terms with different values of k vanish. Integrating over 

transverse coordinates over the iris hole, and assuming Aco/co <?c 1 small 
enough to neglect the dependence on co in Eq. (83) one finds the energy of 
the radiation pulse as 
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e 2 L 2 c , 



ncoa^ 



lk 



xsinc 



L(_co_ 

4 \cy 2 



+ 



coa 2 



i ^ (vi*) 



exp| 



v^oc 3 ^ 2 



(84) 



For the THz wavelength range and for the XFEL spent electron beam ener- 
gies, the constant Lco/(cy 2 ) <sc 1 and can be neglected in all cases of interest. 



6 Radiation of an electron wiggling in an iris line 



Let us now consider the case of an electron emitting radiation from a planar 
undulator with period A w in an iris line. We will be interested in frequencies 
near the fundamental harmonic co r = 2k w cy 2 , where y z = y / Vl + K 2 /2, 
k w = 2n/A w and K is the undulator parameter. We can specify "how near" 
the frequency of interest co is compared to co r by introducing a detuning 
parameter C defined as 



C-— k 

COr 



(85) 



where Aco = co - co r . It can be shown that within the limit 



|C| « k w , 



(86) 



many simplification arise and we can neglect all non-resonant components 
in Maxwell equations. This amounts to neglecting the gradient component 
in the vector g, and to neglecting the constrained motion f*(z) in the Dirac 
S-function, that is, effectively, in the Green's function. Following these pre- 
scriptions, g can be written as 



8 = 



4ne 



exp 



i I dz- 



co 



2y 2 z (z)c 



(87) 



Eq. (51) yields the following expression for the horizontal field component 
E: 



E = 



ecodcA 



s^JJ 



L w /2 a In 

j dz'exp(iCz') j dr' J dcp'G(r,r',(p,(P',z-z')5(r') , (88) 



-W2 
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where O s = K/y, L w is the undulator length and A ;J = / [K 2 /(4 + 2K 2 )] - 
}i[K 2 /(4: + 2K 2 )]. Inspection of Eq. (65) shows that only the term with n = 
yields non-zero contributions. One therefore obtains the following azimuthal- 
symmetric expression for the field at the exit of the undulator: 



E(r) = - 



/(l - 2A)ed s A n 



oo 



icvUl-2A)T 



Lw/2 



X 



r licviAl-2A)z' s 
I tfe'exp(iCz')exp \ 2 



2coa 2 2 ) 
(89) 



Performing the integration in dz' finally yields 



E(r) = - 



i(l-2A)ed s Ajj 
ca 2 



oo . 



icv 2 (l 



2A) L, 



2o«? 2 



xL 7 „sinc 



2coa 2 




(90) 



Eq. (90) can be used as a basis to calculate the energy per unit spectral interval 
per unit surface at the sample position. In our case of interest, the distance 
between sample and extracting mirror, L s , is much longer compared with the 
length of the undulator. Assuming L s » L w and accounting for coherence 
through the form factor F(co), the field energy per unit spectral interval per 
unit surface at the sample position is given by 



dW 
dcodS 



N 2 e \F(co)f 



4tt 2 

e 2 Lld 2 A 2 



,2^ 1 [ v 2 p c 3/2 b 1/2 



xsmc 



2a 2 co) 



(91) 



It should be noted that in the calculation of the square modulus of the 

field, |e| , crossed terms with different values of k vanish. Integrating over 
transverse coordinates over the iris hole, and assuming Aco/co ~ 1/N W <sc 1 
small enough to neglect the dependence on co in Eq. (91) one finds the energy 
per spectral interval Aco/co as 



e 2 L 2 ..co6 2 A 2 , l2 ~ i 



( v^W 2 
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xsinc 




Aco 

CO 



(92) 



As expected, different modes do not interfere. The total pulse energy is given 
by the sum of the energy of each of the empty iris guide eigenmode excited 
within the open waveguide. 

One figure of merit of interest is the ratio between the pulse energy gener- 
ated and transported within the iris guide, and the pulse energy generated 
and transported in free space. In the space-frequency domain, the field dis- 
tribution from the undulator in free space is easily derived. One should 
take the limit for a — > oo in Eq. (88). Accounting for Eq. (73), we get the 
following well-known expression in the far zone limit (i.e. for z » L w ): 



~ , ^ 6sOJeL w A }} 
x (Z ' ) = 2c^ 6Xp 



.co6 2 z 



2c 



sine 




e x , 



(93) 



where = r/z defines the observation direction, and 6 = \6\. The distribution 
of energy radiated per unit solid angle, per unit frequency interval is defined 
as 



dW 



cz^ 



■N 2 e \F(co)f 



dcodD. 47i 2 

Integration over angles can be easily performed leading to 



(94) 



e 2 L w co 2 6 2 A 2 , l7 Am 
W„= 8c2 ' " Nl\F(<o)\ 2 ^. (95, 

This expression corresponds to the particular choice of detuning parameter 
C = 0, at which the pulse energy (in free space) achieves its maximum. 
The ratio between the energy of the pulse generated and transmitted with 
quasi-optical techniques and that of the pulse generated and transmitted 
in free space provides a measure of the departure from ideal performance. 
Accounting for Eq. (32), Eq. (92), and Eq (95) we get the following expression 
for this ratio: 



W 2 ch- u 



W fe 



7i coa z 



xsmc 



cv 



exp 



co^ 



T C+ 2P^ 



Ok \ 



(96) 
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THz radiation source at the LCLS baseline 
phase I THz edge radiation source 



radiation energy losses per one iris 0.1 I 
total energy losses in 100 m - long 
transport line 25 1 



15 m -long vacuum 
chamber with iris line 


THz radiation 
to the access 
maze 






electron beam 
to the beam 
dump 




■ i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 


















Beryllium diameter of the vacuum chamber 14 cm Beryllium 

foil distance between iris 30 cm foil 

diameter of the hole 1 1 cm 



Fig. 23. Scheme of generation of edge radiation at the LCLS baseline. 



THz edge radiation source 



exit of baseline 
undulator 



34 m 



Thermal 

Barrier THz radiation 
682 m to access maze 



15 m -long iris line 



vacuum chamber 



foil 



electron beam 
to beam dump 



foil 



Fig. 24. Installation of the iris line for the THz edge radiation source downstream 
of the LCLS baseline undulator. 

7 Scheme for generating THz edge radiation from the LCLS baseline 



The THz edge radiation source proposed in this paper is compatible with 
the layout of the LCLS baseline and can be realized with minimal additional 
efforts. The vacuum chamber equipped with iris line and outcoupling sys- 
tem can be installed in the unoccupied straight vacuum line behind the 
baseline undulator, Fig. 23 and Fig. 24. The transport of the THz beam in 
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THz radiation 




Fig. 25. Scheme of transport of the THz radiation to the experimental X- ray floor. 

double electron bunch generation for THz pump - X-ray probe experiments at the LCLS baseline 



peak current 

a 100 ns delay 




1 nC 0.25 nC 



Fig. 26. Double electron bunch generation for THz pump/X-ray probe experiments 
at the LCLS. The THz radiation transport through the access maze introduces a 
path delay with respect to the X-ray path, which has to be compensated with the 
introduction of a delay between two bunches. 

the LCLS baseline case only be done with the use of the dump area access 
maze, which is reproduced schematically in Fig. 25. 

The THz transmission line transporting the THz radiation through the ac- 
cess maze introduces a path-delay of about 30 m. Since THz pump/X-ray 
probe experiments should be enabled, we propose to exploit a two-bunch 



44 



at b >2R 

it is possible to match the 
incident and outgoing 
radiation 



iris line 



incident path 



2a I 
2R 



-| i*— b — »| f 



outgoing path 



ilane mirror 



all mirrors deflect the beam by 90 c 



Fig. 27. Geometry of the transmission line turn. The losses are calculated at the 
wavelength A = 0.1 mm. 

mode of operation in order to cope with the delay between THz and X-ray 
pulses, Fig. 26. The two bunches should be precisely separated by a cer- 
tain temporal delay and such delay should be matched to that introduced 
by the THz transmission line. The main idea is, first, to produce two elec- 
tron bunches at a fixed distance with respect to each other. Second, the 
bunches are accelerated and compressed. The beam energy and position is 
kept nearly identical, while they progress one after the other through the 
undulator and THz setup. A way to create two pulses delayed with respect 
to each other is to split the UV laser pulse which illuminates the cathode, de- 
lay one of the pulses, and then recombine the pulses on the same trajectory 
heading towards the rf gun. One method for accomplishing this is described 
in [22, 24]. A multibunch mode of operation was also used at the SLAC linac 
[25]. Routine operation with about 60 ns spacing was achieved. It should be 
noted that for the present scheme (in contrast with that discussed in [22] ) 
we do not require an essentially flat RF amplitude between bunches. 

The operation of the proposed THz source is insensitive to the emittance 
and energy spread of the electron beam. An analysis of the parameters of 
the THz source shows that it will operate reliably even for emittance and 
energy spread exceeding the LCLS commissioning results by an order of 
magnitude. 

The THz transmission line for the LCLS baseline includes at least six 90 
degrees turns, and will exploit plane mirrors as functional components. If 
the pipe of the transmission line has a diameter smaller than the distance 
between the irises it is possible to match incident and outgoing radiation 
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k=1 mode 
k=2 mode 




2,0x10 12 4,0x10 12 6,0x10 12 8,0x10 12 1,0x10 13 

v[Hz] 



Fig. 28. Edge radiation pulse energy emitted within spectral bandwidth 
Aco/co = 10% and transported at the sample position as a function of frequency. 
Partial contributions of individual modes of the circular iris guide are illustrated. 
The sample is set 100 meters away from the extracting mirror. The curves are cal- 
culated with analytical formulas in Eq. (32) and Eq. (84). Here N e — 6.4 • 10 9 (1 nC), 
L = 15 m, b — 30 cm, 2a — 11 cm. The bunch form factor used is shown in Fig. 7. 

without extra losses in these irregularities, Fig. 27. 

It is possible to calculate the energy in the pulse, assuming a 10% spectral 
bandwidth. The bunch form factor considered here is given in Fig. 7. We 
also consider a total length of the iris line of 100 m. The energy loss in the 
line are accounted for. The energy per pulse at the exit of the iris line as a 
function of the frequency v = co/(2n) is shown in Fig. 28. 

The maximal value of the pulse energy is achieved at A — 0.06 mm. When 
the bandpass filter is tuned to this value of A, the expression for the total 
edge radiation pulse energy at the sample can be written in the form: 



W[mJ] * 0.34 . (97) 

CO 

The energy as a function of frequency exhibits a low frequency cutoff due 
to losses in the transport line and a high frequency cutoff due to form factor 
suppression. 
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THz radiation source at the LCLS baseline 



phase I] THz undulator source 



radiation energy losses per one iris 0.04 1 
total energy losses in transport line 11 1 



vacuum chamber 
with iris line 



undulator period 150 cm 
number of periods 10 



THz radiation 
to the access 
maze 



electron beam 
to the beam 
dump 



i i i i i i i 



i i i i i i i i i i 



i i i i i 



Beryllium 



diameter of the vacuum chamber 14 cm 
distance between iris 30 cm 
diameter of the hole 1 1 cm 



wiggling 
amplitude 
6 mm 



Fig. 29. Scheme for generating THz undulator radiation at the LCLS baseline. The 
installation of the THz undulator does not perturb the edge radiation setup and 
does not interfere with the THz edge radiation source operation. The losses are 
calculated at the wavelength A = 0.1 mm. 



exit of baseline 
undulator 



34 m 



Thermal 
Barrier 
682 m 



THz undulator 
with iris line 

-I 



1 



electron 
beam 



planar geometry 
length of period 1.5 m 
number of periods 10 
Max. field strength 15 kG 



Fig. 30. Installation of the THz undulator downstream of the baseline undulator. 

8 Scheme for generating THz undulator radiation from the LCLS base- 
line 

As discussed in Section 2, an optimal expansion strategy for the LCLS THz 
source should be based on the installation of a THz undulator, Fig. 29, 
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Fig. 30. The undulator proposed is an inexpensive planar electromagnetic 
device with 10 periods, each 1.5 m long. At the electron beam energy around 
5 GeV and at a THz wavelength around 0.2 mm, the peak value of the 
magnetic field is about 1.5 T. Undulator radiation will be produced around 
the fundamental harmonic, and the fundamental harmonic will be tuned by 
adjusting the magnetic field strength. The wavelength range between 0.05 
mm and 0.2 mm provided by this powerful radiation source will overlap 
with a large part of the THz gap. This will allow applications for pump- 
probe experiments combining THz and X-ray radiation 10 . 

Note that in the case of THz undulator source at LCLS the electron beam 
emittance is much smaller than the radiation wavelength. This means that 
the electron beam can safely be modeled as filament beam. Computer codes 
like SRW [26] and SPECTRA [27] can be used to study the characteristics 
of the output radiation in the space-frequency domain up to a wavelength 
where the influence of the vacuum chamber is negligible. However, in the 
case of the THz undulator source at the LCLS, wavelengths of order of 0.1 
mm and 15 m - long undulator yield a radiation diffraction size of order of 10 
cm. This rough estimate indicates that vacuum chamber effects are expected 
to play an important role. In this case, conventional computer codes fail to 
predict the correct radiation characteristics. 

Summing up, in view of practical applications to the THz undulator at 
LCLS, there is a need for a theory of undulator radiation in the presence of 
a waveguide. In [28] we developed a theory of undulator radiation in the 
presence of the vacuum chamber with circularly symmetric cross-section. 
In order to efficiently couple radiation into the iris transmission line, for the 
LCLS case we proposed to use a THz undulator source together with an iris 
line, Fig. 29 and Fig. 30. The task that one has to solve differs from the cir- 
cular vacuum chamber case [28] in the formulation of boundary conditions. 
Using Leontovich boundary conditions gives a good approximation in the 
case of a metallic vacuum pipe with resistive walls, and drastically simpli- 
fies the solution of the electrodynamical problem [28]. The problem of mode 
excitation in an iris waveguide can be solved following the same approach, 
where Leontovich boundary conditions are now substituted by Vainstein 
boundary conditions. In our practical case of interest the wiggling ampli- 
tude of the electrons in the undulator is small with respect to the dimension 
of the iris waveguide, Fig. 30. This greatly simplifies analytical calculations. 
In free-space and under the resonance approximation undulator radiation 
is horizontally polarized. This is actually a replica of the undulator polar- 



10 There is a trade-off between the far-infrared and the X-ray frequency range 
achievable in the case a THz undulator source is used. For instance, at an electron 
energy around 10 GeV the longest far-infrared and shortest X-ray wavelengths 
available for operations are around 0.05 mm and 0.25 nm respectively. 
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k=1 mode 
k=2 mode 
k=3 mode 




0,0 2,0x10 12 4,0x10 12 6,0x10 12 8,0x10 12 1,0x10 13 

v[Hz] 



Fig. 31. Undulator radiation pulse energy emitted with a spectral bandwidth 
Aco/co = 10% and transported at the sample position as a function of frequency. 
Partial contributions of individual modes of the circular iris guide are illustrated. 
The sample is set 100 meters away from the extracting mirror. The curves are calcu- 
lated with the analytical formulas in Eq. (32) and Eq. (92). Here N e = 6.4 • 10 9 (1 nC), 
N w = 10, A w — 1.5 m, b - 30 cm, 2a — 11 cm. The bunch form factor used is shown 
in Fig. 7. Curves correspond to the particular choice of the detuning parameter C 
at which the pulse energy in the second mode achieves its maximum. 

ization properties. Moreover, the field exhibits azimuthal symmetry. This 
properties are the same when an iris line is introduced 11 . 

It is possible to calculate the energy in the pulse, assuming a 10% spectral 
bandwidth. The bunch form factor considered here is given in Fig. 7. We 
also consider a total length of the iris line of 100 m. The energy loss in the 
line are also accounted for. The energy per pulse at the exit of the iris line as 
a function of the frequency v = co/(2n) is shown in Fig. 31. 

It can be easily shown that the energy per spectral interval, Eq. (92), does 
not depend on the energy of the electron beam nor, in our case of interest 
with K » 1, on the undulator parameter. In other words, we assume that 
the fundamental harmonic will be tuned by adjusting the magnetic field 
strength at some fixed electron beam energy, but there is no need to specify 



These properties are lost in the case of metallic pipe when Leontovich boundary 
conditions are introduced [28]. 
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these parameters. 



It is instructive to calculate the ratio W/Wf s in Eq. (96) and find the departure 
of the results presented in Fig. 31 from the free space, ideal performance. 
Calculating the Eq. (96) at the sample position and at the optimal wavelength 
A 0.06 mm, we obtain W/W is - 0.8 . This comparison demonstrates the 
effectiveness of the proposed setup. 

Finally, we note that according to our calculations, at a wavelength around 
0.06 mm, and for ten cycles in the pulse, one obtains a 0.3 GW peak power 
level on the sample (see Fig. 31). With strong focusing, these THz beams 
will approach the high field limit of 1 V/ atomic size. 



9 Conclusions 

In this paper we presented a design of a THz edge radiation source and 
a possible extension to a THz undulator radiation source for the LCLS 
baseline. We began our considerations from the generation of THz radia- 
tion from the spent electron beam downstream of the baseline undulator. 
Transmission of the THz beam can only be accomplished with quasi-optical 
techniques. In this article we proposed to use an open beam waveguide 
such as an iris guide, and we presented a complete theory of iris guides. In 
particular, eigenmodes were studied both numerically and analytically. In 
order to efficiently couple radiation into the iris transmission line, it is desir- 
able to match the spatial pattern of the source radiation to the propagating 
mode of the transmission line. To solve the matching problem, we proposed 
to generate the THz radiation within the iris guide. 

Following a derivation of the iris guide Green's function we presented the 
the first exhaustive theory of edge and undulator radiation within an iris 
guide, and we exemplified it in the case of the proposed THz source at LCLS. 
We specialized our consideration to the case of a circular iris waveguide. 
The electric field was found as a superposition of the iris waveguide modes, 
and was studied for the LCLS parameters case. 

The concepts presented here are applicable not only to the LCLS, but also to 
other XFEL facilities. In fact, we regard the concepts presented here, based 
on the combination of edge radiation and undulator radiation with an open 
waveguide setup, as the start of a novel direction towards the construction of 
THz sources at XFELs. Therefore, the present work should not be considered 
as comprehensive of all possible solutions. On the contrary, this work is 
presented with the goal of stimulating interest and of opening the door to 
new possibilities for THz radiation production and transport at XFELs. 
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